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Abstract 

Conventional cluster and virial expansions are generalized to momentum dependent inter- 
particle potentials. The model with Lorentz contracted hard core potentials is considered, e.g. 
as hadron gas model. A Van der Waals-type model with a temperature dependent excluded 
volume is derived. Lorentz contraction effects at given temperature are stronger for light 
particles and make their effective excluded volume smaller than that of heavy ones. 
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The Van der Waals (VdW) excluded volume model has been used to describe hadron 
yields in relativistic nucleus-nucleus collisions (see e.g. JI], Q and references therein). This 
model treats the hadrons as hard-core spheres and, therefore, takes into account the hadron 
repulsion at short distances. In a relativistic situation one should, however, include the Lorentz 
contraction of the hard core-hadrons. This problem was discussed in the literature (see e.g. 
Ref. 0, 3). In this paper the cluster and virial expansions are generalized to velocity dependent 
inter-particle potentials. This extension is used to construct the VdW model for Lorentz 
contracted rigid spheres which may be used to simulate hadrons. 

The canonical partition function for the gas of classical (Boltzmann) particles takes the 
form 
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where V and T are the system volume and temperature, g is the number of internal degrees 
of freedom (degeneracy factor) of the particles, tu, = {m? + k^)^-^^ is the dispersion relation of 
free particles with masses m. The particle interactions described by the function U in Eq. (|l]) 
are given by the sum over pair potentials: 
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In contrast to the usual statistical mechanic treatment of the pair potentials, the Uij are 
assumed to be both coordinate and momentum dependent Uij = M(rj,kj; Vj.'kj). This gen- 
eralization is necessary, if Lorentz contraction effects of hard spheres are to be taken into 
account. Introducing the Mayer functions 
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Eq. (|T]) can be presented as 
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with the short notation rfxj = gdrid]<.i/{27T)^. Similarly to the standard procedure one can 
introduce the cluster integrals ^ 



UJi 

T 



h = y J rfxi exp j^- 
h = 7^ J rfxidx2 exp 



gT^ (m 



2\V 
1 



12 



(ixi(ix2(ix3 exp 



T 



+ /12/23 + /13/23 + /12/23/13) , 



(/12/: 



13 



(5) 
(6) 
(7) 



and present the canonical partition function in the familiar form 
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where the summation in Eq. (|^) is taken over all sets of non- negative integer numbers {rrii} 
satisfying the condition 
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Note, however, that the cluster integrals defined above are different from those used in standard 
statistical mechanics as here nontrivial momentum integrations are included. Condition 
(P) makes the calculation of (i) rather complicated. This problem can be avoided in the 
grand canonical ensemble: the grand canonical partition function can be calculated explicitly 
[z = exp(/x/T)): 
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In the thermodynamical limit the pressure p and particle number density n are calculated in 
the grand canonical ensemble in terms of the asymptotic values of the cluster integrals: 
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The virial expansion represents the pressure in terms of a series of particle number density 
and takes the form 

oo 
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Substituting p ([TT| ) and n (12) into Eq. (13) and equating the coefficients of each power of z, 
one finds the virial coefficients ai in terms of the cluster integrals 

1 ^2 462 263 .... 

ai = 1 , aa = - ^ , «3 = - , - (14) 

Let us recall, first, the derivation of the standard VdW excluded volume model. Then it 
is extended by adding the Lorentz contraction of the moving particles. Keeping the first two 
terms of the virial expansion ([13D the following result is obtained: 

p{T,n) = Tn (1 + a2n) . (15) 

It is valid for small particle densities (i.e. n << 1/02)- The usual (momentum independent) 
hard core potential for spherical particles with radius Tq is Uij = M(|rj — |). Here the function 
u{r) equals to for r > 2ro and cxd for r < 2ro. The second cluster integral (^ can easily be 
calculated in this case: 

h = - ^T) ^-frl . (16) 

Therefore 02 = 4t>o, where Vo = 4:7irl/3 is the particle hard core volume. The VdW excluded 
volume model is obtained as the extrapolation of Eq. ([T5|) to large particle densities in the 
form 

p(T,n) = . (17) 
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For practical use the pressure is given as a function of T and /i independent variables, i.e. 
in the grand canonical ensemble. This is done by substituting n = {dp/dfi)T into Eq. (p!7D, 
which then turns into a partial differential equation for the function p{T,fi). For the VdW 
model ([T7| ) the solution of this partial differential equation can be presented in the form of a 
transcendental equation 

p(T,/i) = T</)(T)e'^/^exp(-^) = p^T^fi-a^p). (18) 

Eq. (|I^) was first obtained in Ref. |^ using the Laplas transform technique. With p{T, fi) 
(the solution of Eq. (|18|)) the particle number density, entropy density and energy density are 
calculated as {u = fj, — a2p(T, /i), 02 = 4wo): 
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Here the superscripts id in the thermodynamical functions (p!8|-pl|) indicate those of the ideal 
gas. 

The excluded volume effect accounts for the blocked volume of two spheres when they 
touch each other. If hard-sphere particles move with relativistic velocities it is necessary to 
include their Lorentz contraction in the rest frame of the fluid. The model suggested in Ref. 
0] is not satisfactory: the parameter 02 = 4 f o of the VdW excluded volume model is confused 
there with the proper volume of an individual particle - the contraction effect is introduced 
for the proper volume of each particle. In order to get the correct result it is necessary to 
account for the excluded volume of two Lorentz contracted spheres. 

Let Ti and rj be the coordinates of the i-th and j-th particle, respectively, and kj and 
kj be their momenta, fjj denotes the unit vector fj^ = rjj/|rjj| [vij = |rj — rj|). Then for 
a given set of vectors (fjj,kj,kj) for the Lorentz contracted rigid spheres of radius there 
exists the minimum distance between their centers rjj(f j^; kj, kj) = min|rjj|. The dependence 
of the potentials Uij on the coordinates rj, and momenta kj, k^ can be given in terms of the 
minimal distance as follows 
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u{Ti,kf,rj,kj) = i (22) 
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The general approach to the cluster- and virial expansions described above is valid for this 
momentum dependent potential, and it leads to Eqs. (p!7| , p!8| ) with 
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The new feature is the temperature dependence of the excluded volume 02 ( |25| ) which is due 
to the Lorentz contraction of the rigid spheres. The pressure and particle number density are 



still given by Eqs. (pr8|Jl9| ), but with temperature dependent 02 (T) (^). However, Eqs. (|20| , 
2TD are now modified, e.g. 
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In contrast to Eq. ( pTD the energy density (p^ ) contains the extra term which appears also in 
the entropy density. The excluded volume 02 (T) (p3D is always smaller than Avo- It has been 
proven rigorously that a2{T) is a monotonously decreasing function of T and, therefore, the 
additional term in Eq. ( plf ) is always positive. Let us introduce the notation 

a2(T) = 4v,f{T) . (25) 

The function /(T) depends on the T/m ratio. It can be calculated numerically and its behavior 
is shown in Fig. 1. The simple analytical formula 
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is found to be valid with an accuracy of a few percents for all temperatures. The asymptotic 
behavior of /(T) is the following: 1 — 0(T/m) at T << m and c + 0(m/T) at T >> m. 

If one assumes that all types of hadrons have at rest the same hard core radius then the 
Lorentz contraction effect leads to different VdW excluded volumes for moving particles with 
different masses: for light particles (e.g. pions) the excluded volume (at given T) is smaller 
than that for heavy ones. Fig. 1 shows that at T = 150 MeV the value of 02 in the nucleon 
gas (m = 939 MeV) decreases by 10% in comparison to its nonrelativistic value 4fo, whereas 
for pions (m = 140 MeV) 02 shrinks at the same T by almost a factor 2. This is simply 
because light particles are more relativistic than heavy ones at given temperature typical for 
high energy nuclear collisions, T = 120 -i- 170 MeV. 

As an example. Fig. 2 shows the particle number density of the pion gas {fi = 0, g = 3) 
with To = 0.5 fm. The particle number density is calculated according to Eq. ([T9|) for three 
different models: the ideal pion gas (02 = 0), the VdW model with constant excluded volume 
(02 = 4wo) and the VdW model with Lorentz contraction (02 (T) is given by Eq. (^31) ). It can 
be seen from Fig. 2 that at low T the pion density is small and excluded volume corrections are 
unimportant. Therefore, all three models are similar. The situation changes with increasing 
T: the suppression due to the excluded volume effects are large and different for 02 = 4fo and 
02 (T) (^). The ratios of particle number densities and energy densities of the pion gas for 
two versions of the VdW model (02 = 4fo and a2{T) (|23|) ) are shown in Fig. 3 as functions 
of the temperature. From Fig. 3 one can observe the deviations between these two models. 
These deviations increase with temperature. They are larger for the energy density due to the 
additional positive term in Eq. (P^. 

In conclusion, the traditional cluster and virial expansions can be consistently generalized 
to momentum dependent pair potentials. Hard-core potentials with Lorentz contraction effects 
lead to a VdW model with a temperature dependent excluded volume 02 (T) (^). For light 
particles the effect of Lorentz contraction is, evidently, stronger than for heavy ones. Note 
that smaller values of the pion hard core radius r^^ were introduced in Refs. M, M within the 



standard VdW excluded volume model to fit hadron yield data better. The smaller value of 
the pion excluded volume appears as a consequence of stronger Lorentz contraction for light 
particles. 
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Fig. 1. /(T) as the function of the temperature-to-mass ratio. For heavy particles (e.g., 
nucleons m » T) the volume reduction is just a few per cents, whereas for pions (m ~ T) it 
is about 50%. 
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Fig. 2. The particle number density for three models of the pion gas {/j, — 0, g — 3): the 
solid line corresponds to the VdW model of the Lorentz contracted spheres {tq — 0.5 fm), the 

dashed one corresponds to the ideal gas of point-like particles, and the dotted one corresponds 
to the VdW model without Lorentz contraction for the spheres of a constant radius 0.5 fm. 
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Fig. 3. The dashed hne shows the ratio of the particle number densities of the pion gas 
{g = 3, II = 0,ro = 0.5 fm): the VdW model with Lorentz contraction divided by the VdW 
model without Lorentz contraction. The solid line shows a similar ratio for the energy densities. 



